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Abstract 

I extend the two gluon exchange picture of elastic scattering, known as the 
Low-Nussinov or subtractive quark model, to predict cross sections for double 
pomeron exchange processes. In particular, I calculate pp ppqq where the 
qq partons will appear as jets separated from the final p and p by large gaps in 
rapidity. The predicted cross section is large enough that this process should be 
observable at the Fermilab Tevatron and at the Large Hadron Collider. It can 
be distinguished from the background of ordinary jet production by an absence 
of particles produced in the gap regions. 



1 Introduction 



The exchange of a gluon is the simplest interaction between two hadrons in QCD. It 
corresponds, via s-channel unitarity, to an elastic amplitude dominated by exchange 
of two gluons in a state with vacuum quantum numbers — in particular, a color 
singlet. This provides a simple model, known as the Low-Nussinov model for the 
pomeron that governs diffractive scattering at high energy. The energy dependence 
of the model (s^ in the amplitude) is close to the observed behavior (~ Q), so 

the picture is qualitatively reasonable. 

In this paper, we extend the Low-Nussinov model to predict cross sections for 
double pomeron exchange (DPE) processes |^, ^ ^, which are characterized by 
two large rapidity gaps |lOl- It has been suggested P that these processes will 
be observable at the Fermilab Tevatron, and it is important to try to predict their 
cross sections. DPE will also be an interesting subject for study at the Large Hadron 
Collider (Vi= 14TeV). 

We will focus on qq jet production in pp scattering, where the final state contains 
only the two jets and a p and p that carry > 95% of their original momenta. The 
final p and p have transverse momenta < 1 GeV, putting them too close to the beam 
directions (pseudo-rapidity \ri\ > 7) to be seen with present detectors. The absence 
of particles produced outside the two jets (in Lego variables r] and 0), except in the 
region between them because of soft QCD radiation, contrasts strikingly to ordinary 
events — especially those with a hard scattering — so the signature of DPE-produced 
jets will be unmistakable. Meanwhile, the hard scattering amplitude is under control 
in perturbative QCD, so no new parameters are added to the Low-Nussinov picture 
to make the prediction. 

To calculate two gluon exchange, we need a model for the internal color structure 
of the hadrons to which the two gluon system couples. In this paper, the scattering 
hadrons are taken to be qq bound states of effective "quarks" that have spin zero and 
couple to the hadrons by a point coupling. One might of course prefer qqq for the 
wave function of a baryon, and additional qq pairs and gluons are certainly present 
in wave functions for small momentum transfer. Our simple model may nevertheless 
be adequate, since only the distribution of color in the hadron is significant for the 
calculation, and that can be modelled correctly by adjusting parameters to fit elastic 
scattering. Indeed, only the distribution of color as a function of impact parameter 
really matters since the spin 1 gluon "sees" quarks equally, regardless of their lon- 
gitudinal momentum. This justifies the simplicity of using spin quarks. We will 
also try an exponentially falling model for the wave function, which is more realistic. 
The model dependence will be assessed by comparing results obtained using the two 
different wave function types, with a range of choices for their parameters. 

Higher-order effects such as interaction between the exchanged gluons must be 
important at some level, and is evidenced by the deviation from constant total cross 
sections as a function of energy. Interactions could even build a rather conventional 
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Regge trajectory for the pomeron, with physical glueball states on it at positive t, as 
Landshoff has emphasized recently |Q . More-than-two gluon exchange contributions 
are also not negligible, as can be estimated by eikonalizing the two gluon amplitude 
H|. (The contribution beginning with four gluon exchange has sometimes also been 



called double pomeron exchange |]TT|. It must not be confused with the definition 
of DPE used here.) We will neglect these effects both for elastic scattering and jet 
production. Fitting the model parameters to elastic scattering should reduce the 
consequences of this approximation. 

In spite of its simplicity, it is worthwhile to see what this model has to say about jet 
production in double pomeron exchange, which has not until now been calculated — 
even though more exotic DPE processes of heavy quark ^ and Higgs [|12| production 
have been. The calculation will be presented in considerable detail to make clear 
how it could be applied to other double- diffract ive processes. The model could also 
be applied to hard scattering with the exchange of just one soft pomeron, i.e., hard 
scattering in diffractive dissociation. Examples for study would be 7p ^ qqp, which 
can be observed at HERA and has been calculated in somewhat different models 
0' M' 1^1 i ^^'^ single-diffractive production of which has been looked for by 



CDF |T6|. 



An additional theoretical motivation for this work is that the two gluon exchange 
picture provides an explicit model for a "direct" , "coherent" , or "lossless" contribution 

|T^, in which the full energy of the pomeron is available for hard scattering. From 
a theoretical standpoint, such contributions are interesting because they violate the 
QCD factorization rules that have been established for inclusive processes, and that 
are often assumed without proof for the diffractive subset of final states |T^, pn| . 
They appear as an effective "super-hard" term oc 6{x — 1) in the phenomenological 
parton distribution of the pomeron. Suggestive experimental evidence for a coherent 



contribution has been seen by UA8 |21 



2 Elastic Scattering 

The pomeron is believed to arise from diffractive physics, i.e., to be an s-channel uni- 
tarity phenomenon. We therefore want the imaginary part of the two gluon exchange 
amplitude, which can be calculated for elastic scattering from the discontinuity illus- 
trated in Fig. 1: 



2z(27r)4 J (F - m2) [{k + pi - p^f - m^] 

X D^^{pi, k -*p3,k + pi- ps) Df,^{p2, -k Pi, -k+p2- Pi) ■ (1) 

Note that the imaginary part is conveniently found by cutting the diagram into two 
pieces through the possible physical intermediate states. Alternatively, one could 
calculate the amplitude as a Feynman diagram, in which case there would be an 
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additional diagram in which the two gluons cross each other. The real part would 
cancel between these diagrams because the amplitude has even signature and energy- 
dependence cx: s^. 

Eq. (|I|) contains a factor 8 from the sum over gluon colors. A finite mass rrig is 
included in the gluon propagators to suppress contributions from long distance, as 
an approximation to color confinement. (An alternative modification of the gluon 
propagator is discussed in Ref. p^ .) 

Our model for the discontinuity of the gluon-hadron amplitude, illustrated in 
Fig. 2, is 

Df,u{pi, k^p3,k + pi- Pa) = 



X 



X 



2 J (27r)4 
(pi -qi + 52)/. 



m 



{pi -q2 + qi) 



ipi ~ qiY ~ {pi ~ q2y — 

jps -qi + q2)y _ jps -q2 + qi)u 
(P3 - gi)^ - m'^ {P3 - q2y - 



(2) 



where q2 = k + pi — ps, and g and G are couplings of the scalar quark to a gluon and 
to the hadron. An overall factor | from color is included, although it is not actually 
significant because the coupling strength g'^ is taken as a free parameter of the 
model. 

In view of the delta functions, Eq. appears to be a two-dimensional integral. 
In the high energy limit, however, one of the delta functions can be reserved to 
apply to the d'^k integral in Eq. leaving a three-dimensional integral. To see 
this, introduce the light-cone coordinates p± = (po ^Pz)/V2 and work in a Lorentz 
frame such as the center of mass, where pi^ = p^^ and p2~ = Pa,- are large with 
■s = ij>i + P2Y — 2pi+P2- • Then d^q = d'^q±dq+dq- . Introduce Xi = qi+/pi+ and 



X2 = q2+/Pi+ and use 6{q^ — m?) 
The other delta function becomes 

- 



(5(2gi+gi_ 



qL 



m 



to do the gi_ integral. 



6{k.) ^ 5{p,-k) 
2 q2+ X2 



in the high energy limit and we obtain 

D^,u{pi,k p^,k + pi - p^) = (pi)^ iPi)u Sipi-k) T{k±, {ps - pi)_ 



(3) 



(4) 



where 

T{k±, [n-Pi] 



dxi I dx2 X1X2 5{1 — Xi — X2) d^qi 



TC^ JO JO 
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1 



1 
1 



(5) 



Here q2± = {pi + k — qi)± and rfi 



~2 



m 



XiX2M^ with M the mass of the hadron 



and m the mass of the "quark". We have ignored the differences between (pi)^ 
{P3)u and (pa)^ (pa),,, which are non-leading in s. 
The transverse momentum integrals can be carried out to obtain 



(6) 



where 



F{k^,A_ 



47r2 



dx 



x[l — x) 



(,q± + k±- xA_ 



|2 + 



X (1 — x) 



9^ 
2n 



dxxil 



X] 



AB 



log 



xil — x) 



B + 



(7) 



B 



with A = {k± — xA±) ^ and B = yl + Am? /A . The remaining integral over quark 
momentum fraction x can be done numerically by Gauss-Legendre integration. 

A more realistic model can be made by replacing the two energy denominator 
factors of the form [X — M^]"^ in Eq. (|^ by exponentials oc e~^"^. This leads to 



F{k^,A^) = N / rfa;e-^t(^^^-"^^)'+^™'l/2x(i-x) 



(9) 



where the three constants /5, m, and the normalization parametrize the wave 
function in place of the point-coupling model parameters. The exponential form 
mimics the fragility of actual hadrons at low momentum transfer, which is displayed 
by the approximately exponential fall-off of elastic amplitudes at small —t. The 
structure of Eq. @) can be seen better by writing it as 



F{k±, Aa_) = N f dx e-(/^/2) [(^'+4"^')/^ + (ki+4m^)/(l-x) - Al] 



(10) 



where k^ = fc^ and k2 = (A_|_ — k±)'^ are the squared momentum transfers carried by 
each gluon. A saddle-point approximation to Eq. (|10|) is convenient for speeding up 
numerical calculations: 



F(fcx,Ax) 



A^ 



I27rab 



P 



[a + by 



,-(/3/2)[(a+fe)2-A2] 



:iii 
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where a = ykl + Am? and h = yk2+ Am? . 
Returning to Eqs. (^-(1), we have 



I s 



d?k^ [T{k±, A 



with = —(^pi — p3)2 = —t the momentum transfer. The k± integral can be done 
numerically. 

We have three models for the proton wave function, given by Eqs. (0) and 
(pTp, with the latter two equivalent except in computational convenience. On physical 
grounds, we estimate m = 0.5 (or 0.3) and = 0.3 (or 0.14 or 1.0) in GeV = 1 
units, with the alternatives representing an estimate of the possible range. 

The model amplitude is proportional to as expected for spin 1 exchange. It thus 
describes an energy- independent total cross section and elastic da/dt. To choose the 
remaining two parameters in each model, we fit the elastic and total cross sections to 



o"tot = 65 mb and Cci/o-tot = 0.207, which are based on a fit ||23[ to pp elastic scattering 



data p4| at ^/s = 546 GeV. The parameter values are shown in Table I. 

The parameters of the point-coupling model (Eq. (^) are such that 2m is very 
close to M. This can be understood using the uncertainty principle: the large spatial 
extent of the proton, which is responsible for the large elastic slope, is reproduced 
in the model by the possibility of nearly on-shell dissociation of the proton into its 
constituents. A choice like m = M = mg = 0.3, suggested in preliminary work by 
Berera |]25| for example, would instead make cxei/atot = 0.87, which is much too large. 
Equivalently, but independent of the normalization, it would make the average elastic 
slope cr^ot/lGTTCTei equal to 3.8 GeV~^ instead of 16.0 GeV~^. 

The exponentially falling wave function models (Eqs. (p!0|) , (0)) reproduce the 
actual shape of da/dt much better than the point coupling model. These two models 
are very similar, since the second is just a saddle-point approximation to the first, 
with parameters chosen to give the same atot and aei- Our smallest choice mg = 0.14 
begins to have a slope at t = that is too steep, since the slope diverges in the limit 
mg 0. 

The hadronic discontinuity modelled by Fig. 2 is actually a function of 3 scalar 
variables, say fcf^, k2j_, and {ki± + k2±Y. Tuning the model parameters to fit the 
t-dependence of elastic scattering may therefore not be sufficient to determine it. 
(The ambiguity could be reduced by also fitting the electromagnetic form factor 0.) 
However, we will find that the predictions are not extremely sensitive to the form of 
the model. 

Some important aspects of Eq. (|T^ are especially clear in the forward direction, 
where it reduces to 

Id^kA^^^^^^] (13) 
^T^^J \ kl + m^g ' ^ ' 
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The factor F(0, 0) — F{k±,0) represents the response of the hadronic wave function 
to gluon momentum transfers of k± and —k± . F{0, 0) comes from the "diagonal" 
diagrams in Fig. 2, where both gluons hit the same quark hne so there is no k± 
dependence from the wave function. F{k±, 0) comes from the "off-diagonal" diagrams 
in which the gluons hit different quarks. 

Eq. ([T3|) has two different momentum scales: one associated with the gluon prop- 
agator and one associated with the hadronic wave function. The overall dependence 
on k± is set by the fit to elastic scattering, but the relative contributions are not well 
determined. 

First consider the parameters in Table I for Eq. (|T0]) with the smallest assumed 
mg = 0.14 = mjr, which corresponds to relatively long-range color confinement. We 
find F(0, 0) — F{kj_, 0) ^ F(0, 0) (1— e^^^'^^ ). With these parameters, the cancellation 
between diagonal and off-diagonal terms is extremely important: omitting the off- 
diagonal term would increase Utot by a factor > 5 . Meanwhile, the gluon mass is 
rather unimportant: even setting it to zero would increase cxtot by only a factor of 
1.6 . To emphasize the importance of cancellation between contributions in which the 
gluons couple to the same or to different quarks, which follows from color-neutrality 
of the complete hadron, this picture has been called the suhtractive quark model [Q, 
in contrast to the pre-QCD additive quark point of view. 

Now consider instead the parameters in Table I for Eq. (p!0[) with the largest 
assumed value = 1.0 GeV, which corresponds to very short-range confinement. 
We find F(0,0) - F(A;x,0) ^ F(0, 0) (1 - e'^^^^). With this choice of parameters, 
the off-diagonal term is quite unimportant: omitting it would increase atot by only 
10%. This point of view in which the diagonal terms are dominant corresponds to 
the additive quark picture advocated by Donnachie and Landshoff [§, |26| . 

We will find that the DPE predictions are somewhat different for the different 
choices of m^,, so DPE measurements might eventually be used to decide the correct 
point of view. 

3 DPE Production of qq Jets 

Fig. 3 shows a natural extension of the Low-Nussinov model to describe qq jet pro- 
duction in DPE. These diagrams can be expected to dominate all other contributions 
because the hadronic discontinuities D^^ contain phase- coherent sums over physical 
states, like their counterparts in elastic scattering. As in the case of elastic scattering, 
a diagrammatic calculation would include crossed graphs that individually generate 
real parts which cancel in the sum. The discontinuity method is simpler as well as 
being more intuitively related to s-channel unitarity. 
The absorptive part of the amplitude is 
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X [D/s^ipi, k p3, ki) Dp^{p2, -k p4, k2) + 
D^,f3iPl, -ki p3, -k) D^p{p2, -k2 Pi, k)] 
lf,l-{Pb-ki)-iu , lul■{Ph-k2)l^, 



X m(p5 



(P5 - A;2)^ 



(14) 



where ki = pi — p^ + k and = P2~ Pi — k. Note that we use true spin 1/2 massless 
quarks here, in contrast to the effective "quarks" used in the wave function modeL 
Eq. reduces this to an integral over transverse momentum, and the contributions 
from the two sets of diagrams in Fig. 3 are equal in view of the symmetry T{k^, A±) = 
T{A± — k±, A±). This leads to 



M 



^9 

16 vr^ 



d^k^f{k^)A. 



T{k_L,p3±) T{-k_L,pi±) 
{kl + m2) (fc2^ + m2) + m2) 



A, 



<?<? 



U{P5) 



l-pil-{,P5-ki)-i-p2 ^ TP21-{P5 



k2) 1-pi 



{P5 - kiY 



{P5 



k,V 



(15) 
(16) 

V{P6) (17) 



where we set pi± = p2± = for incoming particles in the ± z direction. 



To compute , first compute J2 ^Iqik') Aqq{k) where the sum is over q and 



q helicities and k' is an integration variable independent from k. Neglect non-leading 
powers in s by dropping pi_, p2+, /c2+ . Let 



— * 

Q± = {P5± -P6±)/2 



where is essentially the transverse momentum of each jet since p5± + Pg± = 
~P3i. ~ Pii-i with |p3_|_| and \pai] limited to < 1 GeV by the proton wave function. 
The transverse momentum in the loop integration is also limited by the proton wave 
function, so Q± is large compared to all other transverse momenta. Keeping only the 
leading power in gives 



40 



3QI cosh^5 



i{k') a{k) cosh^ 5 + h{k') h{k) sinh^ 5 ] (19) 



where 



5 = (^5 - y6)/2 



(20) 
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and 



a{k) 
b{k) 



{ki)x {k2)y + {ki)y {k2)x 
{ki)x {k2)x - {ki)y {k2)y 



(21) 
(22) 



with k 



P3±; 



k 



2± 



= —k± — p4±, and Q± taken to be in the x direction. 
Eq. (|1^) includes a factor of y from color and a factor of 5 to sum over the quark jet 
flavors d, u, s, c, b. (Double-diffractive top production will be a welcome newcomer 
at LHC.) 

The cross section is 



da 



\M 



qQ\ 



(i2p3_L d'^p4_L d'^Q± dy5 dy^ 



212 ^8^2 



(23) 



where and y^ are the rapidities of the two jets. Both terms in Eq. (|T9|) contain a 
function of k[^ times a function of k±. Their contributions to lA^g^p can therefore be 
computed as absolute squares of integrals over k±: 



da 



a. 



d'^P3j_ d^pi^ d'^Q± dy5 dye 



Qj_ cosh 6 



{ca cosh^ S + Ch sinh^ 5) 



where 
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Cb 



3(27r)i4 

10 

3(27r)i4 



d'^k± f{k±) a{k±) 
d'k^f{k±)b{k^) 



(24) 

(25) 
(26) 



These results can also be obtained by calculating the individual qq helicity amplitudes. 
The cross sections are equal for helicities (+1/2 —1/2) and (—1/2 +1/2), and zero 
for (+1/2 +1/2) and (-1/2 -1/2). 

The cross section is independent of overall energy s . The dependence on jet 
transverse momentum is the usual dimensional Qj^. The dependence on 5 = {y^ — 
ye)/ 2 is such that the two jets are usually separated by <2 units of rapidity. For large 
S, the cross section falls as e~'^^~^^' which is dictated by Regge arguments for spin 1/2 
exchange. Similarly, there is no dependence on the average rapidity y^^e = (l/s + l/6)/2 
of the jet pair because the gluons have spin 1. 

In the special case where both leading particles have zero transverse momentum 
P3_L = p4± = 0, the cross section is found to be zero as a result of the azimuthal 
angle integrations in Eqs. (^3])-(^). This implies strong correlations between the 
transverse momenta of the leading particles. A "Regge factorization" assumption, 
whereby the cross section is a product of factors for emission of pomerons by the fast 
particles times a cross section for two-pomeron scattering, would be incorrect. It also 
implies that it is dangerous to estimate the DPE cross section on the basis of the 
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pure forward direction, as is done in somewhat different models for heavy quark 
and Higgs production. 

Integrating over the transverse momenta of both quasi-elastically scattered p and 
p, since these cannot be observed in current experiments, gives 

J 2 

- {Ca cosh' S + Ch sinh' 5) . (27) 



d^Q± dy5 dye Qi cosh S 

When Ca and are integrated over the azimuthal angles of p^^ and p4±, they are found 
to become equal. Hence Ca = Cb in Eq. ^2li). Results are shown in Table I for our 
various choices of proton wave function and gluon mass. We find Ca = C;, ~ 1.2 x 10~^ 
with an uncertainty up or down of a factor of 2.5 . Somewhat smaller results than 
that are found for the rather extreme choice = 1.0 . 



The large Q± limit in Eqs. (|T9|)-(|22|) was computed with the help of Mathematica 
It is interesting to compare it with production of spin zero quarks, which is 
simple enough to work out by hand as follows. The spin 1/2 factors in Eq. are 
replaced by the likewise gauge-invariant form 

{ki - 2p5)f, {k2 - 2pQ)y {ki - 2pG)f, {k2 - 2p^)^ 

(h^^ + (K^^ + ■ 

Eq. ( p!7D is replaced by 

_ s (1 - al3) 

^SpinO - TTT- w-i , m \'^^) 

(1 + a) (1 + (3) 

where 

a = (fci -P5)i/2P5~P6+ 

P = {h-P6)l/2p5+P6- (30) 

It suffices to approximate a = e'^ and P = e"'"^ =^ (1 + a) (1 + /3) = 4 cosh' 5 in 
the denominator. Terms of order l/Q^ must be kept in the numerator because a/3 is 
close to 1, leading to 

^Spino = cosh' (5 (^^-L'^i-L ^^"^2^ - fci±-fc2±) (31) 

{kix k2x — kiy k2y 

) (32) 



s 



2Ql cosh' 5 



where the final form is for Q± in the x direction. The dependence on ki± and k2± is 
the same as in Eq. (p2|), so the cross section again goes to zero in the double forward 
limit psj_ = p4± = 0. The dependence of the cross section on rapidity is (cosh^)""^, 
which falls as e~''^^~^^' for large separation as required by Regge theory for spin 
exchange. 
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4 Experimental Considerations 



A measurement that could be made with the CDF or D0 detectors at Fermilab {pp 
at ^/s = 1800 GeV) would require two jets in the central region of pseudo-rapidity, 
say |?75|, \riQ\ < 1.5 . Defining the jets using a cone radius of 0.7 would leave regions 
of at least 2.2 < |?7| < 4.2 in the two "end-cap" parts of the detector, to observe the 
absence of produced hadrons that distinguishes DPE from ordinary hard scattering. 

The fraction of longitudinal momentum retained by the forward incident proton 
is Xi = P3+M+ = 1 - iP6++P6+)/pi+ = 1 - (e^^ + e^^)QJ^s. Similarly, X2 = 
1 — {e"^^ + e~'^^)Q±/ ^/s is the fraction of momentum retained by the backward anti- 
proton. Requiring Xi, X2 > 0.95 defines the DPE region as 

\V5\ < 1-5 

l^el < 1-5 

e^^+e''^ < 0.05v^/gx 

+ e-^« < 0.05 V^/Q± . (33) 



e 



-V5 



The predicted DPE cross section for qq jets is calculated by integrating Eq. ( p?]) 
over the region defined by Eq. (|3^) . In doing this, I take as(Q^) = 127r/ (23 InQ^/A^) 
with A = 0.2 GeV and = Ql/A. The result is 



cr 



(3.516 Ca + 0.632 Cb) fib for jets with > 10 GeV/c 

(0.246 Ca + 0.040 Cb) fib for jets with Q± > 20 GeV/c (34) 

(0.022 Ca + 0.002 Cb) lib for jets with > 30 GeV/c . 



Using Ca = Cb and taking = 1.2 x 10 as a typical estimate from Table I gives 



a = 



4.98 nb for jets with > 10 GeV/c 

0.34 nb for jets with Q± > 20 GeV/c (35) 
0.03 nb for jets with Q± > 30 GeV/c . 



These predictions are uncertain by a factor of 2-3 due to the model dependence 
indicated by the spread of values for Ca- The final predicted cross sections will be 
larger because the contribution from gluon jet production is yet to be calculated. 

One might expect the cross section to be reduced by the following "tmin" effect. 
The four-momentum transfer to the leading proton is 

ti = (Pi-Psf = -IpL + (1 -Xi)'mJ]/Xi (36) 

which becomes — p|x in the Xi ^ 1 limit that is assumed in our calculation. To 
correct for this approximation, the predicted cross section should be reduced by a 
factor fti e^*™ where tmm = ^1 + P3± < and B ~ 16 GeV~^ based on elastic 
scattering. A similar factor would be expected for the anti-proton. However, this 
effect is found to be small enough to neglect in the region Xi, X2 > 0.95 . 
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Ordinary hard scattering generates a background to DPE that I estimate using a 
HERWIG QCD Monte Carlo simulation in the manner described in Ref. The 
predicted cross section for 2 jets, each with > 10 GeV, in the DPE region defined 
by Eq. ( |33D is 30 fih. This cross section is nearly 4 orders of magnitude larger than 
the signal. It is also nearly 10~^ of the entire minimum bias cross section, making it 
much too large to permit experiments to trigger on every such event. 

Imposing a rapidity gap condition on one side, by requiring zero particles of trans- 
verse momentum > 0.2 GeV in the range Max(?75,?76) + 0.7 < rj < 4.2, reduces the 
background by a factor 1/700. This makes it small enough to permit triggering on 
all such "single-gap" jet events. 

Imposing a rapidity gap condition on both sides by requiring the regions —4.2 < 
rj < Mm{ri^,riQ) — 0.7 and Max(?75,?76) -|- 0.7 < r] < 4.2 empty of particles with 
p± > 0.2 GeV leads to a HERWIG-predicted background cross section of 1.0 nb. 

Our predicted cross section for quark jets alone is a factor of 5 larger than this 
background, so the DPE signal should show up clearly as an "extra" contribution at 
zero multiplicity in the particle multiplicity distribution for the gap regions of two-jet 
events. As a further test of the model, the DPE region could be tightened to \ri^\, 
{tiqI < 1.0 or < 0.5, which would increase the minimum rapidity gaps from 2.0 to 
2.5 or 3.0 . This would very strongly decrease the background from zero-multiplicity 
fiuctuations of ordinary jet production. Of course, it would be better to extend the 
observed gap regions to larger \ri\, or still better to detect the leading p and p; but 
those options require additions to the detectors. 

5 Conclusion 

We have combined the two gluon exchange model of the pomeron with leading-order 
perturbative QCD for hard scattering to predict cross sections for pp — > ppqq. The 
process shown in Fig. 3 gives the dominant contribution due to phase-coherence of the 
sums over intermediate states represented by the hadronic discontinuities in Fig. 2. 

A similar calculation of pp ppgg is in progress. It is somewhat more complicated 
because many more diagrams make up the appropriate gauge-invariant set. The only 
anticipated difference from the qq result is that large rapidity separations will be 
possible between the jets, as a result of having spin 1 exchange in place of spin 1/2 
between the jets. This could be observable at the LHC, but only by means of detectors 
with a wider coverage in pseudorapidity than those proposed so far. It cannot be 
observed at the Tevatron energy because jets with a large rapidity separation would 
have too large an invariant mass to be produced in the DPE region. 

Our calculation resembles other QCD predictions, in that it contains a long dis- 
tance scale non-perturbative part (the hadronic discontinuity which is related to a 
wave function), and a short distance scale (high-Q^) part that is calculated pertur- 
batively. It differs from other predictions, however, in that the non-perturbative part 
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has been obtained by fitting to low-Q^ elastic scattering rather than to a different 
high-Q^ process; and in that there is no factorization theorem to guarantee success. 

A special feature of this exclusive DPE process is that, unlike rapidity gaps created 
by other color singlet exchanges, there is presumably no additional "survival proba- 
bility" factor needed to account for gaps that are filled in by incidental exchanges of 
color, e.g., due to additional soft gluons exchanged between the incident beam par- 
ticles. This is like elastic scattering itself. Of course, there will be some suppression 
due to the fact that soft particles from the jets can spread widely from the nominal 
jet axes. The effect of such particles can be reduced somewhat by defining gaps as 
an absence of particles above some threshold like p± > 0.2 GeV/c |]10| . 

The experimental signature of our process is two hadronic jets separated by < 2 
units of rapidity, and back-to-back in azimuthal angle. The final p and p are at such 
small angles with respect to the beam directions as to be undetectable in present 
experiments. Installing "Roman Pot" detectors to cover the very small angle region 
would be valuable because it would easily eliminate all backgrounds to DPE, and 
because there are interesting correlations predicted between the transverse momenta 
of the two leading particles relative to each other and relative to the plane of the jets. 
In particular, the predicted cross section vanishes when both leading particles are at 
zero transverse momentum. This strongly contradicts a naive assumption of Regge 
factorization. 

The predicted cross section for the qq process alone is ~ 5 nb. This is large enough 
to be studied easily at the Tevatron and at the eventual Large Hadron Collider at 
CERN. To make the study, it will be necessary to have an experimental trigger for the 
rapidity gap signature on at least one side of the detector. It will also be necessary 
to use sufficiently low luminosity running that the rapidity gaps are not filled in by 
particles from additional pp collisions that occur during the same beam crossing. 

The production of jets discussed here is only one of many possible DPE processes, 
since gg — >■ qq could be replaced by any other hard scattering with a two gluon initial 
state. Some suggestions are given in Ref. 0. A further dramatic possibility would 
be DPE production of a Higgs boson g, |T^. 

Some preliminary work on the subject of jet production in DPE was presented at 
the Fermilab Small-x Workshop p5|, p9[. 
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Figure Captions 

1. Two gluon exchange ( "Low-Nussinov" ) model for the elastic amplitude. The 
dashed line denotes an s-channel discontinuity. 

2. Quark model for the hadronic discontinuity in Fig. 1. 

3. Two gluon exchange model for Double Pomeron Exchange production of qq jets. 
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